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ABSTRACT 

We have used surface plasmon resonance to inves- 
tigate the nucleic acid binding properties of the core 
protein of hepatitis C virus, a disordered protein 
believed to chaperone the genomic RNA. It was pre- 
viously shown that a peptide (peptide E) corres- 
ponding to the association of two basic clusters of 
core enhances the annealing and the dimerization of 
nucleic acid fragments derived from a stem loop 
(SL2) in the 3' untranslated region of the hepatitis 
C virus genome. However, strong aggregation of 
nucleic acids by core or peptide E in the excess of 
the latter precluded the characterization of their 
binding parameters up to now. By careful design 
of surface plasmon resonance experiments, we 
obtained accurate binding parameters for the inter- 
action of peptide E with SL2-derived oligonucleo- 
tides of different lengths and sequences, in form of 
stem-loop, duplex or strand. Peptide E was found to 
bind in a salt dependent manner to all oligonucleo- 
tides assayed. Affinity data identify at least two 
binding modes, of which one is independent of 
sequence/structure, and the other is specific to 
the SL2 stem-loop fold. Stoichiometry data support 
a multi-motif binding model allowing formation 
of higher-order complexes. We propose that the 
modular binding mode demonstrated for structured 
RNA-binding proteins also applies to this disordered 
chaperone and is relevant to its activity. 

INTRODUCTION 

Hepatitis C virus (HCV) is an enveloped, single-stranded, 
positive-sense RNA virus that belongs to the Flaviviridae 
family and causes liver diseases worldwide. Its 9.6 kb 
genome contains a single open reading frame flanked by 



5' and 3' untranslated regions (UTRs) that play essential 
roles in RNA replication and translation. Post-transla- 
tional cleavage of the polyprotein encoded by the HCV 
genome yields 10 viral proteins, among which is the nucleo- 
capsid protein, named core [for review (1-3)]. The mature 
core is a protein of ~21 kDa with multiple functions in the 
viral hfe cycle. Its highly basic and hydrophilic N-terminal 
domain was shown to bind nucleic acids as well as a number 
of host and cellular proteins. Its C-terminal hydrophobic 
domain mediates interactions with lipids and membrane 
proteins [for review (4—7)]. 

This work focuses on the RNA-binding properties of 
the N-terminal part of core as measured by surface 
plasmon resonance (SPR). Core is believed to bind RNA 
and act as a chaperone, a function also displayed by the 
core proteins of other Flaviviridae (8). Indeed, in vitro ex- 
periments have shown that core facilitates structural re- 
arrangements of the 3'-UTR of the HCV genome, leading 
to RNA dimerization (9,10). Core is also able to enhance 
the annealing of complementary HIV- or HCV-derived 
DNA and RNA sequences (9,1 1,12). Although little infor- 
mation is available on the structural and molecular aspects 
of core-RNA interactions in capsid assembly, the in vitro 
analysis of its chaperone properties has identified nucleic 
acid fragments with lengths between 1 6 and 27 nt that are 
re-structured and therefore presumably recognized by 
core. The aim of our study was to investigate whether 
it is possible to identify and quantify distinct RNA- 
core binding modes, potentially related to its different 
functions. For this purpose, SPR measurements were per- 
formed with the molecules used for investigating chaper- 
one activities (10,12) (Table 1). We used a peptide named 
peptide E, corresponding to the association of two of the 
three basic clusters of the N-terminal domain of core (BD2 
and BD3) and shown to support the dimerization and 
annealing activities (10-12) (Figure 1). The RNA frag- 
ments were derived from the sequence of the second 
stem loop (SL2), predicted in one of the possible conform- 
ations of the X-region of the 3'UTR HCV genome (10). 
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Table 1. Oligonucleotides (ONs) used in this study 



Name 


ON sequence 


MW 


Length (nt) 




deoxy-ON 






dDLS 


5'biotin-TCACGGCTAGCTGTGA-3' 


5302 


16 


dPal2 


5'-TGACGCCTAGCTGTGA-3'biotin 


5467 


16 


dPal2 


5'biotin-TGACGCCTAGCTGTGA-3' 


5302 


16 


dRD 


5'biotin-TGCCGCGTAGATGTAC-3' 


5302 


16 


d(TG)4 


5'-TGTGTGTG-3'biotin 


3041 


8 




ribo-ON 






DLS 


5 ' - UC AC GGCUAGCUGUGA - 3'biotin 


5666 


16 


Pal2 


5 ' - UGACGCCUAGCUGUGA - 3'biotin 


5666 


16 


DLSloop 


5 ' - UC AC GGCUUUCUGUGA - 3'biotin 


5604 


16 


cDLS 


5'-UCACAGCUAGCCGUGA-3'biotin 


5650 


16 


Halt-DLS 


5'-AGCUGUGA-3'biotin 


3119 


8 




5'-UCACGGCUAGCUGUGA - 3'biotin 






duplex DLS-DLS 


3'biotin - AGUGUCGAUC GGC ACU - 5 ' 


11332 


16x2 




5'-TCACGGCTAGCTGTGA-3' 






duplex ULcj-qULS 


3'biotin-AGUGUCGAUCGGCACU 


10536 


16x2 




5'-UCACGGCUAGCUGUGA- 3'biotin 






duplex DLS-cDLS 


3'biotin-AGUGCCGAUCGACACU-5' 


11315 


16x2 


SL2 


5' -UCACGGCUAGCUGUGAAAGGUCCGUGA- 3 'bi o t in 


9258 


27 


SL2stem 


5' -AGACGGCUAGCUGUGAAAGGUCCGUCU- 3'bi 0 1 in 


9258 


27 


SL21oop 


5' -UCACGGCUAGCUGUGUUUGGUCCGUGA- 3'bi o t in 


9189 


27 




5' -UCACGGCUAGCUGUGAAAGGUCCGUGA- 3 'bi o t in 






duplex SL2-CSL2 


3'-AGUGCCGAUCGACACUUUCCAGGCACU-5' 


17826 


27 X 2 



The bases that are substituted in the native sequence are marked in bold. The differences in base pairing between the three DLS duplexes are 
underlined. MW, molecular weight; nt, nucleotides. 



Quantitative experimental investigation of pro- 
tein-nucleic acid interactions is challenging because 
specific and non-specific binding events may co-exist and 
be of similar affinity. Although non-specific protein- 
protein interactions are unwanted, non-specific protein- 
nucleic acid binding is Hkely to be an important compo- 
nent of the recognition process, e.g. by faciUtating the 
search for specific sites, a phenomenon amply studied in 
the case of DNA (13-15). Non-specific recognition can 
also be requested for nucleic acid coating in viral 
assembly. A relatively indiscriminate binding mode was 
also suggested to be one of two possible mechanisms for 
chaperoning RNA folding (16). In the first mechanism, the 
protein recognizes and stabilizes a specific RNA conform- 
ation required in the folding pathway. In the second, mis- 
folded intermediates are destabilized by weak non-specific 
interactions, mainly due to electrostatic attraction 
between the positively charged amino acids of the chaper- 
one and the negatively charged phosphates of the RNA 
backbone. To physically separate specific and non-specific 
sites on small RNA fragments in an attempt to isolate 
distinct binding modes is unrealistic. Binding may 
depend on structural features that are not necessarily 
preserved in short fragments. Larger fragments on the 
other hand provide opportunities for various folds and 
for co-existing binding modes that complicate data inter- 
pretation, but better account for biologically relevant rec- 
ognition processes. 

We expect to observe non-specific binding in the inter- 
action system studied here, as peptide E is strongly 



positively charged (pi = 12.5, 18 Arg and 3 Lys in 59 
residues; Figure 1). Furthermore, the N-terminal half of 
HCV core is believed to be disordered, which favors a 
non-discriminatory binding mode (17,18). Disorder, 
often found in chaperone proteins, was proposed to con- 
tribute to recognition diversity by facihtating their adap- 
tation to various nucleic acid sequences and structures 
(19,20). Several activities of core may nevertheless rely 
on the preferential recognition of well-defined RNA 
sequence/structure elements. In particular, the analysis 
of HCV RNA dimerization in the presence of peptide E 
suggested that it is mediated by a 16nt long palindromic 
sequence called dimer linkage sequence (DLS) located 
within SL2, and that the RNA sequence and structure 
influence the core-mediated dimerization (10,21). The pal- 
indromic DLS sequence in its stem-loop and duplex 
forms, together with mutants thereof that are unable to 
dimerize (Figure 1), represent a simple model to investi- 
gate whether peptide E discriminates between ohgonucleo- 
tides (ONs) in stem-loop, duplex and unfolded forms. The 
larger SL2 stem loop used in anneaUng experiments (12) 
was included in our study because it contains structural 
features not present in DLS (Figure 1) that could provide 
other core recognition motifs. 

In addition to challenges inherent to protein-nucleic 
acid binding, peptide E was found to aggregate in the 
presence of ONs, thus preventing binding studies by 
most solution assays. Surface methods represent a 
promising informative alternative because one of the 
partners is fixed so that events subsequent to the binding 
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A Oligonucleotides 



dDLS 




dPal2 




dRD 


DLS 


DLS loop 


cDLS 




Pal2 


T A 




T A 




T A 


U A 


U U 


U A 
C G 




U A 


C G 


C 


G 


G 


G 


C G 


C U 


C 


G 


G-C 


C 


C 


C 


A 


G-C 


G-C 


G C 


C 


C 


G T 


G 


T 


G 


T 


G-U 


G-U 


A C 


G 


U 


C-G 


C 


G 


C 


G 


C-G 


C-G 


C-G 


C 


G 


A-T 


A 


T 


C 


T 


A-U 


A-U 


A-U 


A 


U 


C-G 


G 


G 


G 


A 


C-G 


C-G 


C-G 


G 


G 


T-A 


T 


A 


T 


C 


U-A 


U-A 


U-A 


U 


A 



SL2_a SL2_b SL2 stem SL2loop_a SL2loop_b SL2 loop_c 

G A G A G U 

U A U A u G U U 

GAyGU^ GA UG UU GU 

U-G c A U-G c U C-G U-G 

C-G G A C-G G-U. G-U C-G 

G-U A A G-U A-Jl A-U.. G-U 

C I ^-^u C I C-G C-G c I 

G-C G-C G-C G-C^ G-C^ G-C 

G-C G-C G-C G-C G-C G-C 

C-G C-G C-G C-G C-G C-G 



G- 


C 


G- 


C 


C- 


G 


A- 


U 


G- 


C 


A- 


U 



A-U A-U A-U A-U A-U A-U 

C-G C-G G-C C-G C-G C-G 

U-A U-A A-U U-A U-A U-A 

B Peptides 

HCVCore protein 

m Istnpkpqrk tkrntnrrpq dvkI fpgggqi vnnvYTjipi^^tjP^tBl^^BrWSSSS ^BBBHMSrBtSBrSTI 

BD1 BD2 

SaBB TWAOPn YPWPLYGNEG MGWAGWLLSP RGSRPSWGPT DPRRRSRNLG KVIDTLTCGF ADLMGYIPLV 

BD3 

GAPLGGAARA LAHGVRVLED GVNYATGNLP GCSFSIFLLA L-LSCLTIPAS A 

Peptide D 
GRGSRPSWGP TDPRRRSRNL GK 



Peptide E 



RRGPRLGVRA TRKTSERSQP RGRRQPIPKA RRPEG 



SGRGS RPSWGPTDPR RRSRNLGKG 



Figure 1. Oligonucleotide structures (A) and peptide sequences (B). (A) Secondary structures of the 16mer DLS derivatives (upper row) and the 
27mer SL2 derivatives (lower row) were predicted with the Unafold Web Server program (http://mfold.rna.albany.edu/). The substituted positions in 
the variants are boxed in grey. The 2-3 most stable configurations predicted for SL2 (SL2_a,b) and SL21oop (SL21oop_a,b,c) are shown. (B) Primary 
sequences of the mature HCV core, peptide D (MW 2494 Da) and Peptide E (MW 6754 Da). The three basic domains of the N-terminal region of 
core (BDl, BD2 and BD3) are boxed. 



Step may be prevented. SPR represents a powerful 
approach to investigate protein-nucleic acid recognition 
(22-24). Binding is recorded in real time, and the resulting 
kinetic curves provide information not only on kinetics 
and affinity but also on the binding heterogeneity. 
Furthermore SPR responses are related to the mass 
of molecules at the sensor surface, thus giving access to 
stoichiometry information (25-27). Our data showed 
that peptide E bound, in a salt dependent manner, to 
all ONs assayed, irrespective of their fold and sequence. 
At 150mM NaCl, its affinity (Kq) for most ONs was 
in the range of 0. 1-0.4 |iM. The heterogeneous kinetics 
and Kd around 50 nM observed with SL2 sug- 
gest the existence of a different, sequence/structure 
specific binding mode. Stoichiometrics ranged from 0.7 
to 2.7 peptide E/ON and were related with ON length. 



indicating that the peptide tends to coat the ONs. 
Analysis of the size-stoichiometry correlation led us to 
propose a multi-motif binding model in which two 
peptide E motifs bind two motifs located on the same or 
on different ONs. 



MATERIALS AND METHODS 

CM4 sensor chips and Biacore consumables were obtained 
from GE Healthcare. Streptavidin (ref S4762) was 
purchased from Sigma. 

ONs and core peptides 

The (3' or 5')-biotinylated ONs were purchased from IBA 
GmbH Nucleic Acids product Supply (Gottingen, 
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Germany) and Eurogentec S.A. (Liege, Belgium) and 
purified by the manufacturers by reverse phase high per- 
formance hquid chromatography. Their names and se- 
quences are hsted in Table 1. Peptides E and D were 
synthesized by sohd phase peptide synthesis as described 
in (11). Their sequences are shown in Figure 1. 

Biosensor experiments 

All SPR experiments were performed with a Biacore 
2000® instrument (GE-Healtlicare Biacore, Uppsala, 
Sweden). The running buffer (HEPES-buffered sahne, 
HBS) composition was lOmM Hepes (pH 7.4), 150mM 
NaCl, 3.4 niM EDTA, 0.005% P20 surfactant. (3'- or 5'-) 
labelled biotin-ONs (Table 1) were immobilized on 
streptavidin-coated CM4 sensor chips. CM4 sensor chips 
were preconditioned by duphcate injections of 10 mM 
HCl, 50 mM NaOH, 0.1% (w/v) SDS, each with a 
contact time of 10 s, and water with a contact time of 
20 s. Before covalent immobilization of streptavidin, 
traces of biotinylated products that could remain in the 
flow system were neutralized by injecting a streptavidin 
solution (0.1 mg/ml in HBS) for 5-6 min through all flow 
cells. Streptavidin was then immobilized using standard 
amine-coupHng methods. The surface was activated with 
a mixture of 0.5 mM l-ethyl-3 -(3-dimethylaminopropyl) 
carbodiimide hydrochloride and 0.2 mM N- 
hydroxysuccinimide for 7 min, followed by a 6 min injec- 
tion of streptavidin at 10-100 |ig/ml in lOmM sodium 
acetate (pH 5) and then blocked with 1 M ethylenediamine 
hydrochloride (pH 8.5) (instead of ethanolamine hydro- 
chloride) for 7 min. Finally, the surface was subjected to 
four pulses (100 nl) of 1 M NaCl, 50 mM NaOH at a flow 
rate 50^1/min to wash out all non-covalently bound 
streptavidin. This coupling method resulted in a density 
of 2000-3000 resonance units (RU) of streptavidin on all 
four flow cells. Before immobilization, the single-strand 
ONs, diluted at a low concentration of 10 nM in HBS, 
were heated at 90° C for 1 min and put on ice before injec- 
tion onto the streptavidin-coated chip at 10 |il/min at 25°C 
to generate a 10-30 RU ON surface. ON duplexes were 
preformed by mixing equimolar quantities at high concen- 
trations (>25 |.iM in HBS) of the two ONs (except for the 
duplex DLS-DLS where 50 ^M of DLS were used), 
heating the mix at 90° C for 1 min and then allowing it 
to reach slowly room temperature. In the case of SL2- 
cSL2, DLS-dDLS, the non-biotinylated ON was added 
in 20% excess. A dilution at 10 nM of the duplexes was 
used for injection onto the streptavidin-coated chip to 
reach a density <45 RU. Peptides were serially diluted 
in HBS (containing 0.01% P20 surfactant and NaCl as 
specified) to the final concentrations indicated in Figures 
2 and 3 and injected at a flow rate of 50 (il/min for 1 min. 
The responses were also monitored during a 3-min dissoci- 
ation phase. Although most binding experiments were 
characterized by a fast dissociation, the surfaces were sys- 
tematically treated, after the 3-min dissociation step, with 
0.1% SDS followed by 1 M NaCl in water, each for 30 s, 
to remove any residual protein. Each concentration series 
was run in duplicate (and partially in triplicate or quad- 
ruplicate). Afl binding data were collected at 25°C. 



Data processing and steady-state analysis 

Data processing and analysis was performed using 
Scrubber2 (BioLogic Software Pty Ltd; biologic.com.au). 
The responses were double-referenced (28) by subtracting 
from the signal on the ON surface, both the signal 
recorded on the reference surface (streptavidin without 
ON) and the signal of a blank injection with running 
buffer. The steady-state values were plotted against the 
peptide concentration. The Kq values were determined 
by non-Unear regression analysis of the steady-state 
values using a single site model: [ON. Pep] = Rmax-exp x 
[Pep]/(KD + [Pep]), where Rmax-exp is the maximal 
binding capacity, [ON. Pep], the concentration of 
complex at steady state, [Pep] the concentration of the 
free peptide and K^, the equilibrium dissociation 
constant. 

The binding stoichiometry was calculated as the ratio of 
experimental to theoretical maximal responses (25-27). 
The experimental maximal response, Rmax-exp (SPR 
response achieved when the immobilized binder is 
saturated with the injected binder), was fitted from 
affinity data. The theoretical maximal response, 
Rmax-theor. Can bc Calculated because the SPR signal is 
directly related to the mass of molecules at the sensor 
surface. Rmax-theor was predicted for a one:one binding 
mode (each ON molecule binds one peptide E molecule): 

Rmax-theor = MWpep X RUon/MWqn, whcrC MWpep, 

MWoN are the molecular weights of peptide and ON, 
respectively, and RUqn the SPR signal corresponding to 
the amount of immobilized ON. The binding stoichiom- 
etry was expressed as number of peptide molecules bound 
per ON (Stpep/oN) or as number of nucleotides bound per 

peptide (St„t/Pep)- Stpep/ON = Rmax-exp/R max-theor? Stnt/ 

Pep = total number of nucleotides (in ON)/Stpep/oN- 

The distance between two ONs in the dextran matrix 
was calculated as foUows. The surface density (RU) is 
related to the mass/mm^ by the relation 1 RU = 1 pg/ 
mm^ and was translated into a volume density (mass/ 
vol) assuming a thickness of lOOnm of the dextran layer. 
The number of molecules/vol is calculated assuming a 
homogeneous distribution. This allows to estimate the 
volume of the cube that contains one molecule: 
(MWoN X 10'')/(RUoN X 6) A^ and therefore the distance 
between two ONs [see also (26)]. 

RESULTS 

Deoxy-ON binding properties of peptide E 

We first optimized the conditions of SPR experiments 
using deoxy-ONs instead of ribo-ONs. To investigate 
whether peptide E is able to discriminate between deoxy- 
ONs with different folds, we compared its interaction with 
dDLS and with two variants thereof, in which the palin- 
drome is disrupted: dPal2 and a random sequence called 
dRD. Their sequences and folds are shown in Table 1 and 
Figure 1, respectively. Initially, the peptide was 
immobilized on the sensor surface by amine coupHng 
and flushed with the ON. This assay design should 
prevent potential peptide aggregation and allow the 
analysis of various ONs on a same peptide surface. 
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NaCI ISOmM 



A 1-0 r 



DLS-17RU_Rmax 20 

B 0.50 



NaCI 200 mM 




0 20 40 60 80 100 120 0.10 

Pal2-23RU Rmax27 



1.0 1.5 2.0 



D 15 

1.2 
0.9 



E 0.9 





1.0 1.5 2.0 



Naa_150ml\A 



o 



duplex DLScDLS-29RU_Rmax 1 7 duplex SL2cSL2-44RU_Rmax 1 7 




0 20 40 60 80 100 120 



0 20 40 60 80 100 120 



1.0 1.5 2.0 



NaCI_150mM 



Naa_200mM 



SL2-33RU Rmax24 




0 20 40 60 80 100 120 

Time (s) 



0 20 40 60 

Time (s) 



0.5 1.0 1.5 2.0 

Concentration (|jM) 



Figure 3. Typical sets of SPR binding curves and equilibrium fits for the interaction of peptide E with ribo-ONs. The name of the ONs, together 
with their capture levels (RU) and expected saturation responses (R„,iix-iheor) for 1:1 binding, are indicated above each set of curves. Experiments were 
performed at 150 (A, D, G, H and J) or 200 (B, E and K) mM NaCI. The colour code is as in Figure 2. The concentration series of peptide E 
solutions were injected over the surfaces in duplicate or quadruplicate. The equilibrium fits corresponding to data shown in the left and central panels 
are shown in the right-hand panels (C, F, I and L), and resulting binding parameters are listed in Table 3. 
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However, binding was only observed at a relatively high 
peptide E density (>300RU) and at suboptinial sah con- 
centration for carboxyniethyldextran sensor surfaces 
([NaCl] < 50mM). No affinity parameters could be 
determined owing to complex kinetic profiles and to pro- 
gressive loss of binding activity during repetitive cycling. 
In contrast, when the deoxy-ONs were immobilized (by 
biotin-streptavidin capture) and flushed with peptide E, 
highly reproducible sensorgrams that allowed steady- 
state determination of Kq were obtained at low deoxy- 
ONs density (Rmax<30RU) in the presence of 150mM 
NaCl, after double-referencing of the data. We tested with 
dPal2 that the capture performed by their 3' or 5' ends 
generated similar responses. Typical sets of SPR binding 
curves are shown in Figure 2A-D. The Kq and stoichi- 
ometry (Stpep/oN) values derived from fits of the equilib- 
rium response versus peptide E concentration 
(Figure 2A'-D') are reported in Table 2. The peptide 
was found to bind to the stem-loop dDLS with an equi- 
hbrium affinity (Kd) of 0.2 |iM and a stoichiometry (Stpep/ 
on) of 0.9 peptide/ON, and to the single stranded dPal2 
and dRD with a 4-6-fold weaker affinity (0.9-1.2 |.iM) and 
2-fold higher stoichiometry (2-2.2 peptides/ON). The 
binding parameters measured for dDLS and its variants 
thus indicate that peptide E interacts differently with the 
stem-loop dDLS as compared with the single-strands 
dPal2 and dRD. 

We also recorded binding curves using the shorter 
peptide D that corresponds to the basic domain 3 (BD3) 
of core, or using deoxy-ON fragments unrelated to the 
SL2 sequence. At 150mM NaCl, binding of peptide D 
to dDLS was undetectable (Figure 2E), and the signal 
observed on dPal2 (Figure 2F) and dRD (not shown) 
surfaces was weak. No or httle binding was detected 
when peptide E was injected over surfaces containing the 
shorter deoxy-ONs d(TG)4 (Figure 2G) and d(A)8 (not 



Table 2. Equilibrium parameters of peptide E binding to deoxy-ONs 



Name 


Kd (hM) 


Stoichiometry 


Number 










of 






Peptide/ 


Nucleotides/ 


replicates 






ON 


peptide 








Stpep/ON 


Stnl/Pep 






NaCl 150mM 








dDLS 


0.24 ± 0.05 


0.9 ± 0.1 


17 ± 3 


12 


dPal2" 


0.9 ± 0.3 


2.2 ± 0.2 


7 ± 1 


5 


dRD 


1.2 ± 0.4 


2.0 ± 0.6 


8 ± 2 


4 


d(TG)4 


Not detectable 










NaCl 50 mM 








dDLS'' 


0.04 


1.6 


10 


1 


dRD'' 


0.04 


1.9 


8 


1 


d(TG)4" 


1.1 


2.2 


4 


1 



"The binding parameters cannot be precisely evaluated because the 

highest peptide E concentration used is below 10 x Kd 

''In 25 mM Tris-HCl (pH 8) instead of 10 mM Hepes (pH 7.5). 

Kd values were calculated by non-linear regression analysis of the data 

in Figure 2A'-D' (and Supplementary Figures SI and S2) using the 

Langmuir binding model. The standard deviations were calculated on 

the mean from >4 replicates. The stoichiometry (Stpep/oN) was obtained 

from the ratio Rmax-exp/Rmax-theor. 



shown), which may result from their small size or their 
sequence unrelated to DLS. 

In an attempt to evaluate the electrostatic contribution 
to complex formation, binding experiments were per- 
formed in the NaCl 300 to 30itiM concentration range. 
No binding of peptide E to dDLS could be detected 
at 300 mM NaCl (Figure 2H). At 100 mM NaCl, the 
post-injection phases indicated a slightly slower dissoci- 
ation (not shown), a phenomenon ainplified at 50 mM 
NaCl (Figure 21 and J, Table 2). Binding curves at the 
lowest salt concentrations follow coinplex kinetics. 
Nevertheless, equihbrium values indicate a higher 
complex stabihty with decreasing [NaCl] (Table 2). 
Stoichiometry increases with decreasing salt in the case 
of dDLS but not dRD so that differences in binding 
properties between the stem-loop dDLS and the single- 
strand dRD becoine undetectable. Furthermore, binding 
was also detected at low salt concentration when peptide E 
was injected over d(TG)4 (Figure 2K and L: 50 and 
30 mM NaCl, respectively) and d(A)8 surfaces, or when 
peptide D was injected over a dDLS surface (not shown). 
Thus, the salt-dependence of binding indicates a strong 
electrostatic contribution, as expected for the interaction 
between ONs and a positively charged peptide. 

Binding of peptide E to DLS and its variants 

Having optimized experimental conditions using deoxy- 
ONs, we analysed the binding of peptide E to the corres- 
ponding ribo-ONs. Because differences in peptide E 
binding to deoxy-ONs were most visible at 150mM 
NaCl as compared with [NaCl] < 150 mM, and because 
binding was undetectable at 300 mM NaCl, the experi- 
ments with ribo-ONs were performed at 150 and 
200 niM NaCl. 

Eight DLS-derived ribo-ONs were assayed (Table 1): 
the wild-type DLS, Pal2 (carrying two replacements in 
the stem), DLSloop (carrying two replacements in the 
loop), the complementary sequence cDLS and three 
duplexes: DLS with itself (duplex DLS-DLS), with 
dDLS (duplex DLS-dDLS) and with cDLS (duplex 
DLS-cDLS). An 8mer corresponding to half of the DLS 
sequence (half-DLS) was also assayed. 

Typical binding curves recorded with DLS and Pal2 at 
150 and 200 mM NaCl are shown in Figure 3 A, B, D and 
E, and the corresponding fits in Figure 3C and F (see also 
Supplementary Figures SI and S2). Binding parameters 
are reported in Table 3. Binding curves typical of DLS- 
derived duplexes are shown in Figure 3G. Peptide E was 
observed to bind all of these 16mer ribo-ONs with similar 
affinities (Kq = 0.2-0.4 ^M). This differs from what has 
been observed with the deoxy-ONs: peptide E bound 
dDLS with a 4-6-fold stronger affinity compared with 
dPal2 and dRD (Table 2). A lower affinity was 
measured for half-DLS (Kd = 3.1hM). Stoichiometrics 
were variable with a trend similar to that observed with 
16mer deoxy-ONs: less peptide was bound to the stem- 
loops (Stpep/oN ~1 for DLS and dDLS), as compared 
with the strands Pal2, dPal2 and dRD (Stpep/NA between 
1.5 and 2). The variant stem-loops DLSloop and cDLS 
showed stoichiometrics of '^1.5 and 0.7, respectively. The 
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Table 3. Equilibrium parameters of peptide E binding to ribo-ONs 



Name 


Kd (hM) 




Stoichiometry 


Number of 










replicates 






Peptide/ON 


Nucleotides/peptide 








Stpep/ON 


Stnl/pep 






NaCl ISOmM 










u.zu ^ u.uo 


1 .U It u.z 


1 / i J 


1 n 

1 u 


Pal2 


0 1 0 + n 08 

U. 1 7 ^ U.UO 


1 5 ± 0 1 


1 1 ± 1 




DLSloop 


0.44 ± 0.10 


1.4 ± 0.1 


11 ± 1 


5 


cDLS 


0.25 ± 0.02 


0.7 ± 0.03 


24 ± 1 


3 


Half DLS" 


3.1 ± 0.1 


0.9 ± 0.2 


10 ± 2 


4 






1.6 ± 0.04 


21 ± 0.5 


2 


diinlpx DT S-dDI S 


0 17 ± 0 05 


1.5 ± 0.05 


21 ± 0.7 


2 


duplex DLS-cDLS 


0.29 ± 0.05 


1.4 ± 0.01 


23 ± 0.1 


2 


diinlex SI ''-cSI 7 


0 33 ± 0 05 


2.7 ± 0.3 


20 ± 2 


2 


SL2 


0 05 ± 0 01 


1.8 ± 0.2 


16 ± 2 


9 




0 05 + 0 01 


1.7 ± 0.1 


16 ± 1 


2 


SI 7lnon 


0 05 ± 0 01 


1.6 ± 0.1 


17 ± 1 


2 




NaCl 200 mM 








DLS" 


2.1 ± 0.2 


0.9 ± 0.02 


18 ± 1 


2 


Pal2^ 


1.5 ± 1.1 


1.2 ± 0.2 


13 ± 2 


3 


DLSloop" 


2.7 


1.0 


17 


1 


cDLS" 


4.9 


0.8 


20 


1 


duplex SL2-CSL2" 


4.3 


2.2 


25 


1 


SL2 


0.24 ± 0.19 


1.3 ± 0.1 


20 ± 1 


7 


SL2stem 


0.34 ± 0.05 


1.3 ± 0.01 


20 ± 0.1 


2 


SL21oop 


0.38 ± 0.14 


1.2 ± 0.1 


23 ± 1 


2 



"The binding parameters cannot be precisely evaluated because the highest peptide E concentration used is below 10 x Kd. 

Ko values were calculated by a non linear regression analysis of the steady-state data (see Figure 3 C, F, I, L and Supplementary Figures SI and S2), 
assuming a Langmuir binding model. The standard deviation was calculated on the mean for the number of replicates indicated in the last column. 



responses were found to be specific to peptide E, as no 
binding signal was observed in the presence of non-related 
basic polypeptides such as K14R (sequence KRTAM 
FQDPQERPR; pi = 10.7), scFvQ34S (29) (pi = 10) or 
lysozyme (pi = 9.2) at similar concentrations (not 
shown). Binding parameters measured at 200 mM NaCl 
are also reported in Table 3. In hne with a strong contri- 
bution of electrostatic interactions, equilibrium affinities 
were found to be about one order of magnitude lower at 
200 mM compared with 150mM NaCl (see also 
Supplementary Figure S2). Kinetic fits are reported in 
Supplementary Figures SI and S2 and the resulting 
binding parameters in Supplementary Table SI. The 
determined kinetically as koff/kon was found to be in good 
agreeinent with the Kd value evaluated from equilibrium 
responses reported in Table 3. Binding stoichiometrics 
were either unaffected or slightly decreased by the 
increase in salt concentration. 

Binding of peptide E to SL2 and its variants 

We then investigated the interaction of peptide E with the 
27mer SL2 and two SL2 variants with replacements in the 
stem (SL2stem) and loop (SL21oop) used in anneahng ex- 
periments (12). SL2 and its variants are predicted to adopt 
various folds, of which the most stable, denoted a, b or c, 
are shown in Figure 1. SL2 was also used in the form of a 
duplex with its complementary strand (duplex SL2-cSL2, 
Table 1). Binding parameters are reported in Table 3, and 
illustrative sets of curves are shown in Figure 3H (for the 



duplex at NaCl 150mM), and J andK (for SL2 at 150 and 
200 inM NaCl, respectively). The corresponding equihb- 
rium fits are shown in Figure 31 and L (and 
Supplementary Figures SI and S2). The interaction 
strength between peptide E and the SL2-cSL2 duplex 
was similar to that of the shorter DLS duplexes 
(0.3 ^M), but with a higher stoichiometry (compare 
Figure 3G and 3H). In contrast, SL2 showed a 4-fold 
stronger equihbrium affinity (Kd = 0.05 ^M) compared 
with that of DLS (Kd = 0.2 ^iM). The kinetics at 
150mM NaCl (Figure 3J) were complex, resembhng 
those observed with dDLS or dRD at 50 mM NaCl 
(Figure 21 and J). An evaluation of the equilibrium data 
using a two-site binding model yielded a high affinity com- 
ponent, with parameters close to those obtained with the 
simple Langmuir model (Kd = 0.03 ^M; Stpep/oN = 1-4) 
and a poorly defined (~|iM) low affinity component (not 
shown). No attempt was made to further analyse the 
binding curves because the conformational heterogeneity 
of SL2, also mentioned by Palau et al. (30), precluded the 
study of structure-binding relations. Affinity and stoichi- 
ometry parameters for peptide E interacting with the 
variants SL2stem and SL21oop (Langmuir binding 
model) were comparable with those observed with SL2. 
The SL2stem was not expected to behave differently 
from SL2. The similar KdS measured with SL2 and 
SL21oop indicates that the loop predicted in SL2 is not 
the sole determinant of strong binding. Binding param- 
eters are influenced by the salt concentration, as 
observed for DLS and its variants. 
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DISCUSSION 

Our aim was to investigate for the first time the ON 
binding properties and determinants of peptide E, a 
peptide encompassing two of the three basic domains of 
the HCV core protein and showing strong nucleic acid 
chaperone properties. The inherent disorder of the 
N-terminal end of HCV core comphcates structural 
studies and its tendency to aggregate in the presence of 
ONs impedes binding studies in solution. Thus, we have 
used a surface-based approach to analyse the interaction 
of peptide E with ONs. Peptide E was shown previously to 
increase the dimerization of DLS (10) and the annealing of 
DLS, SL2 and SL2 variants with the corresponding com- 
plementary sequences (12). If these stem-loop structures 
present features that support specific recognition, their 
binding properties to peptide E should be different from 
those of strands and duplexes. We carefully designed and 
optimized SPR experiments to obtain precise affinity and 
stoichiometry parameters for the interaction of peptide E 
with ONs of different lengths and sequences derived from 
the SL2 sequence. Our data show that peptide E was able 
to recognize all ONs analysed. We searched for relation- 
ships between sequence/fold and the binding parameters 
recorded at ISOniM NaCl, for which the largest number 
of measurements was available. 

Affinity and kinetics of the peptide E-ON interactions 

A histogram representation of Kq values obtained at 
150mM NaCl for ONs classified according to their size 
within a given fold (strand, duplex and stem-loop) is 
shown in Figure 4A. No absolute correlation can be 
found between the Kd values and ON size. The 27mer 
stem-loops showed stronger peptide E binding (Kq 
~50nM) compared with the smaller 16mer stem-loops 
(Kd ~0.2-0.4 |aM), and the 16mer strands stronger 
binding than the shorter 8mer strand (half-DLS), but the 
16 and 27mer duplexes displayed similar peptide E binding 
in the 0.2-0.3 |iM range. Also, the three 16mer strands 
(Pal2, dPal2, dRD) behaved differently, with up to 6- 
fold differences in Kq. A correlation between affinity 
and ON form is also not obvious. The 27mers in stem- 
loop form showed a stronger peptide E affinity (SL2, 
SL21oop, SL2stem: Kq ~50nM) compared with the 
27mer in duplex form (SL2-cSL2: Kd ~0.2 ^M), but all 
16mers, whether in stem-loop (DLS, dDLS, cDLS, 
DLSloop), duplex (DLS-DLS, DLS-dDLS, DLS-cDLS) 
or strand (Pal2) form, showed similar KdS in the 
0.2-0.4 (iM range at 150mM NaCl. 

The stronger equilibrium affinity measured for the inter- 
action of peptide E with the 27mer stem-loops (SL2, 
SL2stem, SL21oop) hkely represents a specific binding 
mode, supported by sequence/structure features not 
present in the SL2-cSL2 duplex or in the other ONs. 
The shape of the kinetic curves for the peptide E - SL2 
interaction (Figure 3J) reflects complexities that prevent 
rehable mathematical fitting. The initial binding rates 
tend to be constant, as typically seen when limited by 
mass transport. Simulations (not shown) indicate that 
substantial mass transport cannot be achieved at the low 
ON densities used, when considering the binding 
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Figure 4. Relationship between the ON structure and its binding par- 
ameters to peptide E. The affinity (A), stoichiometry expressed as Stpep/ 
ON (B) or Stnt/pep (C) are represented as histograms for each group of 
ON form (strands, duplexes, stem-loops) classified in order of 
increasing size. 



parameters (Kd and saturation response, Rmax) observed 
at equilibrium. A possible interpretation of the kinetic 
curves is that peptide E initially recognizes a few nucleo- 
tides (high density of small sites) with fast on-rates, before 
re-organizing into the complexes found at equihbrium that 
involve more nucleotides (low density of larger sites). 
Little selectivity and fast kinetics allow the samphng of a 
large number of configurations, which may be relevant in 
facilitating specific binding or in maximizing nucleic acid 
coating. The binding curves for the peptide E-SL2 inter- 
actions indicate multiple successive and/or co-existing 
events that appear to include a distinct, strong-affinity 
binding mode not observed with other ONs, most likely 
promoted by bulges or loops. 

Most other ONs showed similar KdS in the 0.2-0.4 
range at 150mM NaCl. The higher Kd measured with 
half-DLS, dPal2 and dRD (1-3 |.tM range) could not be 
evaluated with high accuracy because equilibrium re- 
sponses at the highest concentrations used are sometimes 
far away from surface saturation. Nevertheless, their 
peptide E binding affinity is clearly weaker than that 
measured with other ONs. This could be explained by 
the smaller size of half-DLS or by the difference in 
chemical structure between deoxy- and ribo-nucleotides 
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and will be discussed further later in the text. The similar 
KqS measured at 150mM NaCl for ONs with different 
folds, together with the salt dependence of binding and 
stoichiometries above one (see later in the text), suggest 
a non-selective, electrostatically driven, binding mode. It is 
reasonable to assume that positive charges of peptide E 
(pi = 12.5) recognize ON backbone phosphates in a 
sequence and structure-independent manner. 

Stoichiometry of the peptide E-ON interactions 

In contrast to Kq values, the binding stoichiometry (Stpep/ 
on) was related to ON size. A histogram representation of 
Stpep/oN for ONs classified according to their size within a 
given fold (strands, stem-loops and duplexes) is shown in 
Figure 4B. The number of peptide E molecules bound per 
ON increased with their length within each group of ON 
conformers. This observation suggests that peptide E 
occupies a defined surface on the ONs, irrespective of 
their sequence. The number of nucleotides statistically 
occupied (Stnt/pep; column « Stoichiometry: Nucleotides/ 
peptide » in Tables 2 and 3) was calculated by dividing 
the length of the ON by Stpep/oN (i-S- if one peptide binds 
to a 16mer, a maximum of 16nt are needed for binding; if 
two peptides bind to a 16mer, a maximum of 8nt are 
needed for binding). Stnt/pep is an alternative expression 
of stoichiometry that is independent of the ON size. 
Stoichiometries (Stpep/oN. Stnt/pep) are calculated 
assuming that the surface-fixed ONs are all capable of 
binding the peptide. It is well known that the covalent 
coupling of proteins can lead to their partial inactivation 
and that, at high densities, a fraction of the proteins may 
become inaccessible to the injected binding partner. Small 
synthetic molecules that are coupled in an oriented manner 
tend to be fully active, as observed by us for synthetic 
peptide binders (data not shown). The ONs captured on 
streptavidin surfaces through a terminal biotin are likely to 
be fully active and accessible, particularly at the very low 
coupling densities used here. This hypothesis is supported 
by the observation that similar Stnt/pep values were found 
for surfaces with different ON density and for ONs with 
different lengths (e.g. all duplexes). A second assumption is 
that the refractive index increment for proteins and nucleic 
acids is the same, as found by Di Primo and Lebars (31), 
who also summarized the divergent conclusions published 
by others. Whatever the exact relationship between RU 
and mass of nucleic acid at the sensor surface, the 
stoichiometries can be ranked and compared. 

The histogram representation of Stnt/pep (Figure 4C) 
indicates that the number of nucleotides statistically 
occupied is relatively constant within each group, but 
differs across groups: 7-11, 20-23 and 16-17 nt are 
needed for peptide E binding to strands, duplexes and 
stem-loops, respectively. The two exceptions within stem- 
loops are the variants DLSloop and cDLS that differ from 
DLS in loop sequence and size, respectively. These excep- 
tions hint to different mechanisms of complex formation 
with the native and non-native 16mer stem-loops, and to a 
possible sequence/structure-dependent binding, despite 
similar affinities. Within strands, the smaller Stnt/pep 
measured for deoxy- compared with ribo-ON strands 



(7-8 and 1 1 nt, respectively) can tentatively be related to 
the larger distance between consecutive phosphate groups 
in extended DN A (6.9 A) as compared with extended RNA 
(5.5 A). The theoretical extended lengths of a 8-nt deoxy- 
and a 1 1-nt ribo-ON are similar (55 and 60 A, respectively) 
and could correspond to the surface occupied by peptide E. 
The relationship between size and stoichiometry suggests 
that peptide E tends to coat the ONs. 

Model of peptide E interaction with strands and duplexes 

The known overaU structure of strands and duplexes 
aUows speculating on the topology of regions occupied 
by the peptide. The stoichiometry (Stpep/ow) of 1.5 
peptide/ON calculated for Pal2 and DLS-derived 
duplexes (Table 3) could indicate the coexistence of differ- 
ent types of complexes on the sensor surface (i.e. mixture 
of 1:1 and 2:1 complexes), or the presence of a fraction of 
inactive or inaccessible ON. However, such phenomena 
are not expected to be systematically reproducible and 
are thus in variance with the highly reproducible value 
of 1.5 Pep/ON consistently measured for Pal2 (5 measure- 
ments) and for the three 16mer duplexes (1.5, 1.6, 1.4 in 
duphcate measurements). This stoichiometry can be inter- 
preted as the formation of higher-order complexes con- 
taining three peptides and two ONs. The simplest model 
allowing the formation of such complexes requires 2 and 3 
binding motifs on peptide E and the ONs, respectively: 
two motifs in peptide E could then coat two ON motifs, 
located on the same or on different ONs. The size of the 
ON motif covered by each peptide motif can be deduced 
from Stnt/pep. In case of Pal2, Stnt/pep is ~llnt. Each 
peptide motif would in this case cover 5-6 nt. The obser- 
vation that peptide E needs twice as many nucleotides to 
bind 16mer duplexes (Stnt/pep = 20-23 nt) as compared 
with Pal2, suggests that a duplex motif contains 5-6 
base pairs. Various configurations of higher-order 
complexes are possible, of which two are drawn for a 
single strand in Figure 5A (left and central panels). 
Similar configurations are possible with duplexes as 
illustrated in Figure 5A (right panel). 

Binding properties measured for the larger SL2-cSL2 
duplex and for other strands do not contradict the 
model. The 27mer SL2-cSL2 duplex is predicted to carry 
five motifs (each containing 5-6 bp) and to be able to form 
5 Pep/2 ONs complexes (Stpep/oN = 2.5) with a Kq in the 
same range as that measured for 16mer duplexes. 
Experimental data are consistent with this prediction 
(Stpep/oN = 2.7 ± 0.3; Kq = 0.33 ± 0.05 |.iM). The param- 
eters measured for half-DLS, dPal2 and dRD differ from 
those discussed so far by nearly integer stoichiometries 
and weaker affinities. Although they are not as accurate 
as those measured for other ONs, affinities are undoubt- 
edly weaker, suggesting that the peptide motif can recog- 
nize <5-6nt, as illustrated in Figure 5B (left panel). 
Indeed, half-DLS is too short to contain two 5-6-nt 
motifs. Similarly, a stoichiometry of 2 for dPal2 and 
dRD would be consistent with a model in which each 
peptide coats at the maximum 8nt (on average 4nt per 
peptide E motif; Figure 5B, right panel), explaining its 
reduced affinity (0.9-1.2 i^M). Although a multi-motif 
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A Two protein motifs ; ON motifs : 5 to 6 nucleotides or base pairs 

Stpep/ON = 1 .5 

single strand duplex 






B Two protein motifs ; ON motifs : < 5 nucleotides 

Stpep/ON = 1 Stpep/ON = 2 





Figure 5. Model for the peptide E-ON interactions. (A) Two among ttie possible configurations for complexes formed of three peptides and two 
16mer ribo-ON strands (left and central panels), valid also for duplexes (right panel): the peptide (dark grey) and ONs (light grey) contain two and 
three motifs, respectively. (B) Weaker binding may result from the presence of fewer negative charges in the ON motifs because of a smaller interface 
(half-DLS) or because of the more expanded deoxy-ON strand compared to the ribo-ON strand. Bridging (not drawn) cannot be excluded when 
stoichiometrics are integers. 



binding model that allows the bridging of ONs by peptide 
E was proposed based on non-integer Stpep/oN (16mer 
strand and duplexes), an integer Stpep/oN does not 
exclude the possibihty of bridging (not drawn in 
Figure 5B). For example, Stpep/oN = 1 if the two peptide 
motifs bind to the same ON, or if the two motifs in two 
peptides bridge two ONs. The mixed (stem and loop) com- 
position of stem-loops, together with the structural and 
binding heterogeneity of SL2, make it difficult to relate 
structure with binding, but multi-motif binding and 
bridging is hkely to apply in this case also. Altogether, 
our data suggest that with a few exceptions, peptide E 
maximizes ON coating, rather than searching for an 
appropriate sequence for getting the highest affinity. 

SPR has been applied previously to the detailed analysis 
of nucleic acid binding by another nucleic acid chaperone, 
the HIV-1 nucleocapsid characterized by two highly 
folded zinc finger motifs (27,32). The data indicated the 
formation of higher-order complexes with an 8mer dON 
at high ON immobilization levels, whereas the same 
ON-nucleocapsid interaction showed a 1:1 stoichiometry 
at low immobilization levels, suggesting that stoichiometry 
depended on the ON density on the chip. An evaluation of 
the average distance between ONs in the dextran matrix 
indicates that the low ON density used in our experiments 
does not prevent bridging of ONs by the peptide. This 
separation was evaluated to be ~350^00A at immobil- 
ization levels in the range 20^5 RU, used in our study. As 
the dextran chains are not cross-hnked and have an 
average length of 100 nm, this is more than twice the 
required bridging distances of a few hundred Angstroms. 
We also observe different Stpep/oN for different ONs, but 
at similar densities. We conclude that the formation of 
higher-order complexes is possible. 



CONCLUSION 

To our knowledge, this is the first characterization of the 
interaction of a disordered RNA-binding protein with 
ONs by SPR. Affinity data identify different modes 
of ON binding. Its sequence/structure-independent over- 
stoichionietric binding to strands and duplexes points to a 
largely non-specific interaction that may be representative 
of a nucleic acid coating activity. A more specific mode, 
possibly related to core chaperone activities, was observed 
with stem-loops, particularly SL2. Such a co-existence 
of both specific and non-specific binding, associated with 
different roles during the viral hfe cycle has already 
been described for the HIV-1 nucleocapsid chaperone 
(33). Stoichiometry data are consistent with a model in 
which two peptide E motifs recognize two nucleic acid 
motifs, allowing the formation of higher-order complexes 
if these are located on different ONs. The presence of two 
motifs on the peptide is consistent with the fact that it is 
composed of two basic domains of core. The ability of 
peptide E to bridge nucleic acid fragments is unlikely to 
be restricted to strands and duplexes (which served to 
develop the model) and has functional significance for 
nucleic acid chaperone activities such as dimerization or 
anneahng, in bringing nucleic acid stretches together. It 
has been largely shown that structured RNA-binding 
proteins are modular and that this modularity is relevant 
to function [for review, see (34,35)]. The modular binding 
mode demonstrated for structured RNA-binding proteins 
is likely to apply also to the disordered HCV core. 
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Supplementary Data are available at NAR Online. 
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